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Massive bigravity models are interesting alternatives to standard cosmology. In most cases, how¬ 
ever, these models have been studied for a simplified scenario in which both metrics take homo¬ 
geneous and isotropic forms [Friedmann-Lemaitre-Robertson-Walker (FLRW)] with the same spa¬ 
tial curvatures. The interest to consider more general geometries arises, in particular, in view of 
the difficulty so far encountered in building stable cosmological solutions with homogeneous and 
isotropic metrics. Here we consider a number of cases in which the two metrics take more general 
forms, namely FLRW with different spatial curvatures—Lemaitre, Lemaitre-Tolman-Bondi (LTB), 
and Bianchi I—as well as cases where only one metric is linearly perturbed. We discuss possible 
consistent combinations and find that only some special cases of FLRW-Lemaitre, LTB-LTB, and 
FLRW-Bianchi I combinations give consistent, nontrivial solutions. 
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I. INTRODUCTION 

The standard ACDM model of cosmology is based on four main assumptions: general relativity (GR) is the correct 
description of gravitational interactions at energies below the Planck scale, the Universe is homogenous and isotropic 
on large scales (the cosmological principle), the energy content of the Universe is mainly in the form of cold dark 
matter (CDM) and a nondynamical cosmological constant A, and all the structure that we see around us originated 
from nearly Gaussian, adiabatic, and scale-independent quantum fluctuations at early times. All these assumptions 
have been tested with high precision using various cosmological data and seem to be in excellent agreement with all 
existing observations. There are, however, various theoretical reasons why one may want to go beyond this standard 
framework. In particular, the assumption that the late-time acceleration of the Universe is due to a cosmological 
constant term has been strongly questioned from the theoretical point of view, as its small but nonzero value preferred 
by observations cannot be explained by fundamental physics [l]. It is, therefore, important and quite natural to ask 
whether the cosmic acceleration can be explained by a different mechanism than a pure cosmological constant. One 
particular possibility, which has attracted remarkable attention over the last decade, is that a modification of GR on 
very large scales might be responsible for the acceleration (see Refs. Si for comprehensive reviews). One of the 
interesting such infrared modifications is to assume that gravitons are not massless as opposed to what GR tells us. A 
nonzero but sufficiently small graviton mass modifies properties of the gravitational interactions on very large scales 
while leaving them indistinguishable from the predictions of standard gravity on small scales where GR is believed to 
be at work. 

GR is a consistent and nonlinear theory of massless gravity and, therefore, has given the possibility of constructing 
various cosmological models. In order to test the implications of massive gravity for cosmology, one similarly needs a 
nonlinear and consistent theory for massive gravitons. Such a theory was, however, not available for more than 70 years 
after the construction of a linear theory of massive gravity by Fierz and Pauli in 1939 Q. This was mainly because 
any attempts at constructing a nonlinear completion of the Fierz and Pauli theory would face a serious obstacle; the 
theory would suffer from the existence of the so-called Boulware-Deser (BD) ghost degrees of freedom [f|, a property 
which would be fatal to the theory. It was only a few years a go t hat a ghost-free and fully nonlinear formulation of 
massive gravity, and its bimetric extension, was constructed ffii ll~il| (see Ref. fl5| for a recent review). The key step for 
this success was to extend the gravitational sector by at least one new spin-2 tensor field with metric-like properties. 
In addition, in order to avoid the BD ghost, the physical metric of the theory has to interact with the new tensor field 
in a very specific way. In the simplest version of the theory, referred to as the de Rham-Gabadadze-Tolley (dRGT) 
theory of massive gravity, only the physical metric, the one which interacts with the matter sector in the standard 
way, is dynamical, i.e. has an Einstein-Hilbert term in the action, while the second metric, often called “reference” 
metric, does not have dynamics. In this case gravitons posses five degrees of freedom. In the bimetric version of the 
theory, referred to as the Hassan-Rosen theory of bigravity, the reference metric is also given dynamics and, therefore, 
gravitons posses seven degrees of freedom, corresponding to one massless and one massive graviton. 
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The dRGT theory of massive gravity has been shown to suffer from a no-go theorem forbidding flat and closed 
Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological solutions on a flat reference metric [ljj]. In addition, 
dRGT with an open FLRW metric or a nonflat reference metric suffers from the so-called Higuchi instability Ez] 
or other types of instabilities [l8l - |23j . One obvious way to avoid the no-go theorem is to give up on exact FLRW 
solutions, i.e. to consider inhomogeneous and/or anisotropic solutions for the metrics. We, however, know that 
the observable Universe on large scales is very close to being homogenous and isotropic and, therefore, non-FLRW 
solutions, if allowed, must not deviate significantly from the FLRW case on observable scales, and must respect the 
observational bounds on inhomogeneity and anisotropy. Non-FLRW effects should be either of very low amplitudes 
or of scales much larger than our horizon so that they cannot be observed. Interestingly, solutions which satisfy these 
conditions have been shown to exist in dRGT [l6|. Other scenarios with non-FLRW solutions in dRGT, either for 
the physical metric or for the reference metric, can be found in Refs. [2ll [fol Gof (see also Ref. (15] for a thorough 
review of the inhomogeneous and anisotropic solutions in massive gravity). The other possibility to avoid the no-go 
theorem and instability issues in dRGT is to extend the theory. An example for such extensions without adding new 
degrees of freedom is the recently proposed generalized massive gravity theory [3ll ] . Another workaround is through 
theories with extra degrees of freedom or violation of certain symmetries; these include for example quasidilaton [32|, 
varying-mass EH [§3[, nonlocal [33 - [37j . and Lorentz-violating [H, 39] massive gravity. Recently, another solution 
to the no-go theorem has been proposed in Refs. [ 40 I-I 42 I] . It has been suggested that the no-go theorem can be 
overcome if at least some matter couples to a hybrid metric, composed of both the physical and reference metrics and 
constructed in a specific way to keep the theory free of the BD ghost up to a cut-off energy scale which is believed to 
be above the strong coupling scale of the theory. This makes phenomenological studies of the theory possible below 
the cut-off scale (see, however, Ref. [Dj for various complications that the cosmology of this theory would need to 
tackle). The revival of the BD ghost is a general feature of the scenarios where both metrics couple simultaneously 
to matter [4(1 [il], IdiUTs! ] ; these include the simplest case where the two metrics couple to matter minimally |~iol [52[ . 

The Hassan-Rosen theory of massive bigravity on the other hand is immune from the no-go theorem and admits 
usual flat FLRW solutions. The cosmology of bigravity has been extensively studied in the literature at both the 
background and perturbative levels (see, e.g., Refs. I53| - |71| ) and in terms of the lensing and dynamical properties of 
local sources [bl] (see also Refs. [72l Pul for cosmological studies of bigravity where matter couples to both metrics 
through a composite metric). It has particularly been shown that the theory can explain the late-time acceleration 
of the Universe in the absence of an explicit cosmological constant and can, therefore, serve as a viable alternative 
to ACDM (5jJ, 63] at the background level. However, when it comes to perturbations, the Hassan-Rosen theory of 
massive bigravity seems to be suffering from various instabilities. Bigravity models can generally be classified into two 
categories, finite and infinite branches. This classification is based on the fact that the ratio of the two scale factors 
for the reference and physical metrics in FLRW solutions is either increasing (finite branch) or decreasing (infinite 
branch) with time, depending on which combinations of the parameters of the theory are nonvanishing. It has been 
shown that scalar perturbations for all finite-branch models, including a simple single-parameter model called minimal 
bigravity model (MBM) Q, are unstable at early times on small scales [56], .64, 66. :68[, @], although the models are 
viable at the background level. The instabilities do not necessarily rule these models out but make their comparison 
to observations difficult as one can no longer employ linear perturbation theory to study their implications for the 
formation of structure. Linear scalar perturbations for the infinite-branch bigravity (IBB), as identified in Ref. Ifi3 ]. 
are on the other hand stable at all times, making the study of structure formation possible for the model. This 
was done in Refs. [67l |68|. where subhorizon scales were analyzed in the quasistatic limit, various modified gravity 
parameters were calculated, and deviations from GR predictions were presented. It was shown that predictions of the 
model are consistent with existing large-scale structure data and the model can be tested by future experiments. IBB 
was, however, shown later to suffer from two other types of instability; it violates the Higuchi bound MM , which 
can potentially be dangerous as the instabilities might appear at higher-order perturbations for nonlinear structure, 
and tensor perturbations have ghost instabilities at early times IzillzfHziJ. Addressing the problem of instabilities 
in bigravity is currently an active field of research. In principle most of the generalizations of the dRGT theory of 
massive gravity enumerated in the previous paragraph can be applied to bigravity to investigate possible resolutions 
to the instability problems and to construct viable alternatives to ACDM. 

Motivated by these instability problems in the standard scenario of massive bigravity with FLRW metrics, in this 
paper we take the first step in exploring one potential route to resolve the obstacles and construct a viable and stable 
model. It is not clear at this stage whether possible solutions are necessarily in modifications of the structure of the 
theory, i.e. at the level of the action. Such possibilities should definitely be explored, but one should also consider 
cases where the structure of the theory remains intact while other classes of solutions are considered. One of these 
possibilities is the class of solutions with non-FLRW metrics for one or both metrics of the theory. As we mentioned 
earlier, this has been shown to be a promising route in the case of dRGT, and it is, therefore, worth investigating for 
bigravity as well. Before studying the cosmological implications of such cases, one needs to check whether solutions to 
field equations exist and whether they are consistent with basic constraints of the theory, such as Bianchi constraints. 
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This is the objective of the present paper. 

Inhomogeneous and anisotropic solutions in bigravity have been studied in the literature. This includes Bianchi 
cosmologies where both metrics are homogeneous and anisotropic 
metric being FLRW and the reference metric being inhomogeneous 


30, Ha], as well as cosmologies with the physical 
261 ]. In this paper we study various combinations 


of FLRW and non-FLRW metrics in a systematic and more general way, and investigate for which combinations 
consistent solutions to the equations of motion exist. Our aim is to identify such combinations without further 
exploration of their implications for cosmology; we leave this for future work. 

The rest of this paper is organized as follows. In Sec. |TT] we review the Hassan-Rosen theory of singly-coupled 
bigravity and present the field equations and Bianchi constraints. In Sec. Mil we study various combinations of 
FLRW, inhomogeneous, and anisotropic solutions for the metrics of the theory. We start this in Sec. IIII Al with the 
case where both metrics are of the FLRW form but have different spatial curvatures. Since, as we will see, this will 
force the metrics to be of the Lemaitre form, we study those solutions in the same section. We move on in Sec. Ml Bl 
with the combinations where one metric is FLRW and the other one is of the Lemaitre-Tolman-Bondi (LTB) form, 
and in Sec. Ml Cl we investigate the solutions where both metrics are LTB. Combinations of FLRW and anisotropic 
but homogeneous metrics are studied in Sec. Ml El in the context of Bianchi type I solutions. In all these sections 
we discuss the consistency of the solutions when matter sources respect or violate the homogeneity or anisotropy 
assumptions. In Sec. m we go beyond the background solutions and investigate the scenarios where one metric is 
perturbed while the other one is kept unperturbed. Our discussions are based on both scalar and tensor perturbations. 
We discuss our results and conclude in Sec. E 


II. THE THEORY OF MASSIVE BIGRAVITY 


The Hassan-Rosen theory of ghost-free, massive bigravity is characterized by the action [13| | 

^■2 /* _ r 

S = —/ d^x\J — det gR g -/ d A x^J— det fRf 

4 

+ m 4 J rf 4 XA/-detgy^/3 n e n fx/g -1 /) + J d 4 x^J- det gC m {g , $), 


( 1 ) 


n—0 


where M g and Mf are Planck masses and R g and Rf are the Ricci scalars for the metrics g and respectively. 
Here g is the standard, physical metric coupled to matter fields $ through the matter Lagrangian £ m , and is 
the reference metric. The action contains five interaction (mass) terms given in terms of five functions e n . These are 
the elementary symmetric polynomials of the eigenvalues of the matrix •*/g -1 /, where \Jg~ x j\J 3 -1 / = g^ ffj.w The 
forms of these polynomials are presented in, e.g., Ref. [l3j| . The quantities f3 n ( n = 0,1, 2,3,4) are free parameters of 
the theory, and m is the mass parameter. In the following we express masses in units of M g and absorb m 4 into the 
parameters /3 n (m is not an independent parameter of the theory). The action then becomes 


i r _ M i [' _ 

5 = — - / d 4 xy/— det g R g -/ d 4 x\J ~ det fRf 

+ J d 4 x\J — detg ^ /3 n e n (V g~ 4 /) + J d 4 x\J- det gC m (g, $). 


( 2 ) 


By varying the action 0 with respect to g^ one obtains the generalized Einstein equation for the physical metric, 


1 3 

1) PngnxY^v (\/g 1 f'j = Tf, Ji,, 


( 3 ) 


n —0 


where R'j JV is the g-metric Ricci tensor, and the matrices Y( n )(A) are defined as 


Y {0 )(X) = I, 

Y W (X)=X-I[X], 

Y(2){X) = X 2 - X[X\ + \l ([A] 2 - [A 2 ]) , 

Y(3)(A) = A 3 - A 2 [A] + ^A ([A] 2 - [A 2 ]) - ij ([A] 3 - 3[A][A 2 ] + 2[A 3 ]), 
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where X = (^\fg , I is the identity matrix, and [...] is the trace operator. 

By varying the action (J2J with respect to the reference metric f gi/ we obtain 


rU - + jp Y^i-irfa-nUxY^ (x/Fv) = o, (4) 

V1 f n =0 

where R is the /-metric Ricci tensor. Under the rescaling —>• MJ 2 f gv , the Ricci scalar Rf transforms as 

Rf —> MjRf, which results in 

V-det fR f -A MJ 2 \J— det fRf. (5) 

The interaction terms in the action then transform as 

4 4 

'Y^P n e n (vF 1 /) -^^2Pne n (M/Vff -1 /) • (6) 

n—0 n —0 

Since the elementary symmetric polynomials e n (X) are of order X n , the rescaling of /by a constant factor MJ 2 
translates into a redefinition of the coupling constants /3 n —> Mjf3 n , which allows us to assume Mf = lQ 

In addition to the equations of motion for the metrics, there are additional constraints on the dynamics of the 
metrics coming from the Bianchi identities and the assumption that the stress-energy-momentum tensor of the matter 
components is conserved, 


3 

2 V 9 (vFv) = 0 , (7) 

n=0 

where V g is the covariant derivative operator with respect to g giy . Any acceptable bigravity solution must satisfy 
the generalized Einstein equations m and Q, as well as the Bianchi constraint CD- In the rest of this paper, we 
investigate various types of the physical and reference metrics, and identify the ones which are consistent with these 
conditions. 


III. CONSISTENCY OF BACKGROUND SOLUTIONS WITH DIFFERENT METRIC COMBINATIONS 
A. FLRW metrics with different spatial curvatures: The need for a Lemaitre reference metric 

We begin our investigation of bigravity with nonstandard metric forms by considering the solutions for which both 
metrics g and are FLRW with generic spatial curvatures k g and kf, where k g , kf = 0, ±1. The usual background 
analysis of the cosmology of bigravity assumes k g = kf. The reason is partly due to the significant simplification of the 
calculations in this case, and in addition it seems intuitively reasonable to assume that the two metrics respect the same 
symmetries and geometries. More complicated cases, however, cannot be excluded a priori. Let us, therefore, leave 
the choices for k g and kf completely generic and study this case in terms of the consistency of cosmological solutions. 
Before we continue we note that this case has previously been studied in Ref. an d shown to be inconsistent using 
the Bianchi constraint. In the following, however, we use the implications of the Bianchi constraint in this case as 
a tool to systematically construct a particular metric combination, FLRW-Lemaitre, which is consistent. Therefore, 
although our results are in agreement with the findings of Ref. (55|, our approach and objectives are different. 

With the assumptions made above the most general forms for the metrics are 

g^dx^dx" = —dt 2 + a 2 (t) dxg 2 , 
f^dx^dx" = -X 2 (i) dt 2 + b 2 (t) dcd) 2 , 

where 

dr 2 

dxg 2 = - -—- + r 2 ( d6 2 + sin 2 ( 9) dcf 2 ) , 

1 — k g r A 

dr 2 

dxf 2 = --—- + r 2 (dd 2 + sin 2 (6) dcji 2 ) , 

1 - k f r 2 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


1 See, however, Ref. |80L which appeared during the completion of this work, for caveats associated with this rescaling. 
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r, 0 and (j> are spherical coordinates, a and b are the scale factors for g^ and respectively, and X is the lapse for 
f gv . Inserting these metric forms into the Bianchi constraint 0 yields the condition 



/3i (2 + k) + 2/3 2 (1 + 2k) —f 3/3 3 k ( — 


3/3i + 2/?2 (2 + k) —h /3 3 (1 + 2 k) I — 


( 12 ) 


on the /-metric lapse X 7 where k = \ l - 9 , and an overdot denotes a derivative with respect to t. Since the 

V 1 “ k f r 

scale factors a and b depend only on time, it is clear from this expression that X cannot be a function of time only 
and is in general a function of both t and r unless k g = kf. This, therefore, shows that FLRW solutions for the two 
metrics with different spatial curvatures are not allowed, already at the level of the Bianchi constraint. 

Now assuming that X is a function of both t and r, Eq. 0 places a constraint on the /-metric scale factor b, 

(P3 r X' + 2 (1 — k) (/3 2 + fcX)) 6“ + 2 (/?2 tX' + (1 — k) (/3i -I- (3 2 X)) ab + /3i rX'a 2 = 0, (13) 

where a prime denotes a derivative with respect to r. We see from this equation that for general choices of the lapse 
X , the scale factor b should also be a function of both t and r. It can be shown that the metrics 0 and 0 with both 
X and b being functions of both r and t cannot be reformulated in FLRW forms by any coordinate transformations 
(see appendix lAl for a detailed proof). We, therefore, conclude that FLRW metrics with different spatial curvatures 
are not consistent. 

Let us now assume that X and b are both functions of r and t. Using the Bianchi constraint 0, we arrive at the 
expressions 


• /3i (2 + At) + 2(3 2 (1 + 2k) —h 3/?3K 

x=l -“- 

3/?i + 2/3 2 (2 + k) —I - Aj (1 + 2k) 

(/? 3 rW' + 2(l —k)(/3 2 + /? 3 X))& 2 

5 2r ((/?! + 

We can, therefore, see that because of different curvatures of reference and physical metrics, the reference metric takes 
a spherically symmetric and inhomogeneous form, where the lapse and all scale factors are functions of both r and t. 
The most general metric forms corresponding to this case are 

g^dx^dx” = — dt 2 + a 2 (t) dx 2 , (16) 

f^udx 9, dx v = — X 2 (t, r) dt 2 + Y 2 ( t , r) dr 2 + Z 2 (t, r ) r 2 dfl 2 , (17) 

where dfl 2 = dd 2 + sin 2 (d)d(j) 2 . The metric in Eq. 03 has a generic spherically symmetric and inhomogeneous 
form. In the literature [U HH, this type of metric is called Lemaitre metric [83|, and the cosmological model built on 
this metric is called Lemaitre model. In GR, the Lemaitre metric arises when we have inhomogeneous matter sources 
H3; in particular when the pressure and density are functions of both temporal and spatial coordinates. In the 
special case of dust or homogenous pressure the Lemaitre metric reduces to the so-called Lemaitre-Tolman-Bondi 
(LTB) metric, where the lapse does not depend on the spatial coordinates and can be rescaled. 

All our arguments so far for the metrics to take the forms (1161) and (1171) when we assume unequal curvatures were 
based only on the Bianchi constraint 0. We can, however, arrive at the same conclusions by analyzing the Einstein 
equations. Assuming k g ^ kf for the metrics with FLRW forms, we have a nonvanishing ^/(l — k g r 2 )/( 1 — kfr 2 ) factor 
in the Einstein equations for the terms corresponding to the interactions between the two metrics. The interaction 
part in the /-metric Einstein equation plays the role of an inhomogeneous source for /^„, which forces it to take a 
Lemaitre form. For the g metric, there is a coupling to the matter source, and by taking an inhomogeneous matter 
source one can in principle cancel the inhomogeneities coming from the interaction terms; as a result, g^ can maintain 
its homogenous FLRW form. This confirms our finding that the metrics should have the forms (1161) and 03 where 
the matter source is inhomogeneous. 

Let us now derive the explicit forms of the Einstein equations for the metrics m and m As argued above, 
we assume that the stress-energy-momentum tensor of the matter source coupled to the physical metric has an 


2 (farX' + (1 — a«) (/3i + fcX)) ab + /?i rX'a 2 


foX) a + (/? 2 + 0 3 X) b) 


(14) 


(15) 
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inhomogeneous perfect-fluid form, 


T g0 = -P (*, r), 

T gi = V (t, r ), 

T g2 = P (*> r )> 

Tg 3 = P (t, r) , (18) 

where p and p are, respectively, the energy density and pressure for the matter source. For simplicity, here we consider 
only the k g = 0 case. The ^-metric Einstein equations for this case read 


3—T + P — 0o + 0i 


(Y + 2Z) (2 Y + Z)Z 


a z 


a 


02 1 


' 03' 


YZ 2 


— + 2— + p — 0o + 0i (X + 2— 
a z a \ a 


L + 2“ + p = A, + A (x + L±l 

a z a \ a 


02 ( 2 

- 02 


Z 2 


a a‘ 
X(Y + Z) 


03 


xz 2 


a* 
YZ\ 

) 


03' 


XYZ 


(19) 

( 20 ) 
( 21 ) 


The equations of motion for /are too unwieldy to be displayed here. One can find the equations for the most 
general case in Ref. (HI). We see from Eqs. (119[) - (121l) that with an appropriate choice of p and p it is in principle 
possible to find a function a (t) which satisfies the equations of motion. 

In conclusion, our bigravity theory does not allow two FLRW metrics with different spatial curvatures, while it is 
possible to have a combination of FLRW and Lemaitre forms for the metrics g gi/ and /^j,, respectively, if the matter 
source takes an inhomogeneous form. In the opposite case of a Lemaitre form for g^ and an FLRW form for the 
equation of motion for will contain a homogenous Einstein tensor part and an inhomogeneous interaction part. 
Since / M „ does not couple to matter, the inhomogeneities cannot be cancelled and, therefore, this metric combination 
in general does not have consistent solutions. 


B. FLRW—LTB and LTB -FLRW combinations 

The next possibility we wish to explore is the case where one of the metrics is FLRW and the other one is LTB. 
For the case where g^ is FLRW and is LTB (we denote this as FLRW- LTB), the line elements have the forms 

g^dx^dx" = — dt 2 + a 2 {t)dx g 2 , (22) 

f^dx^dx" = —X 2 (t)dt 2 + Y 2 (t, r)dr 2 + Z 2 (t, r)r 2 dfl 2 , (23) 

where again dtt 2 = dO 2 +sin 2 (0)d(f> 2 . Here we assume that g^ has a curvature k g . In addition we follow the standard 
recipe for LTB metrics and assume that the physical metric is coupled to a homogeneous perfect-fluid source. In this 
case the (0,0) and (1,1) components of the g^ equation of motion become 

- 3 ad 2 + a 3 ( 0o + p) + a 2 0 1 ^2 Z + Yyjl-r 2 ^ + a ^-3 k g + 0 2 Z 2 + 20 2 YZ^l - k g r 2 ^j + 0 3 YZ 2 ^Jl - k g r 2 = 0, 

(24) 

a 2 + a 2 (p - 0o - 0iX) - 2a ((ft + 0 2 X) Z - a) + k g - 0 2 Z 2 - foXZ 2 = 0. (25) 

Here we have assumed the stress-energy-momentum tensor for the isotropic and homogeneous perfect fluid to be of 
the standard form 


Tgu = (p + p) u%u 0 v + pSI, (26) 

where p = p(t) is the rest energy density of the fluid, p = pit) is its pressure, and u is its isotropic four-velocity. 
It is clear from Eq. (1251) that in general Z cannot be a function of r since all the other quantities in the equation, 
including the /-metric lapse X, are functions only of t. If Z is a function of t only, in order to satisfy Eq. (® Y 
should be of the form Y(t,r ) = A(t)/y/l — k g r 2 , where A(t) is an arbitrary function of t. This then implies that the 
reference metric should also be of an FLRW type with the same curvature k g . 
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In the opposite case, where the physical metric is LTB and the reference metric is FRLW (we denote this as 


LTB -FLRW), we have 

g^dx^dx" = —dt 2 + Y 2 (t,r)dr 2 + Z 2 (t,r)r 2 dfl 2 , (27) 

f^dx^dx 1 ' = —X 2 (t)dt 2 + b 2 (t)dxf 2 . (28) 

The (0,0) and (1,1) components of the equation of motion then read 

X 2 Y (fhZ 2 + 2 p 2 Zb + fob 2 ) yj l~k f r 2 + b (/ 3 2 X 2 Z 2 + 2 p 3 X 2 Zb + p 4 X 2 b 2 - 3 k f X 2 - 3= 0, (29) 

P 4 b 2 X 3 -k f X 3 + (px + p 2 X) X 2 Z 2 + p 3 X 2 b 2 - Xb 2 + 2 Xbb + 2 bX 2 (p 2 + p 3 X) Z - 2 Xbb = 0. (30) 


Similarly to the case for g^, here we again see from Eq. (13(71) that Z in general cannot depend on r. Equation (12U|| , 
therefore, implies that Y should have the form Y(t,r) = B(t)/W 1 — fc/r 2 , where B (t) is again an arbitrary function 
of t. This form for Y then brings the physical metric into an FLRW type, again with the same curvature as the 
reference metric. 

It is important to note that in both metric combinations discussed above the Bianchi constraint is not satisfied, and 
since we get the Bianchi constraint using the covariant conservation of the stress-energy-momentum tensor, we can 
therefore state that in general neither FLRW-LTB nor LTB-FLRW can occur for any choices of covariantly conserved 
T9 

fll/ ’ 


C. LTB—LTB combination 


Let us now study the case where both metrics g^ and v are of LTB forms (we denote this as LTB- LTB), 

g xivdx^ dx u = —dt 2 + A 2 (t,r)dr 2 + B 2 {t,r)dfl 2 , (31) 

f^dx^dx v = -X 2 {t)dt 2 + Y 2 {t , r)dr 2 + Z 2 (t, r)dfl 2 . (32) 

In this case the 0 component of the Bianchi constraint 0 enforces the /-metric lapse X (t) to satisfy the equation 


where 


2 UZ + VY 
2UB + VA ’ 


(33) 


U (t, r) = B(t , r) (/3i A(t, r) + /3 2 Y(t, r)) + Z(t , r) {p 2 A(t, r) + p 3 Y (; t , r)), (34) 

V(t,r) = Px B 2 (t,r) + 2 p 2 B(t,r)Z(t,r) +P 3 Z 2 (t,r). (35) 

The constraint (1351) holds when 

B(t,r) = B(t)s(r), (36) 

Z(t,r) = Z(t)s{r), (37) 

A(t,r) = A(t)q(r), (38) 

Y(t,r) = Y(t)q(r), (39) 

where s and q are functions of r only, obtained through solving Einstein equations. In addition, the (1, 2) component 
of the equation of motion for g^ becomes 

B'A = AB’. (40) 


For fpu the corresponding equation is 

Z'Y = YZ’. 

From these equations we obtain the following relations between B(t), A(t), Z(t ), and Y(t): 


(41) 


B{t)A(t) = A{t)B{t), (42) 

Z(t)Y(t) = Y(t)Z{t). (43) 





We, therefore, find that A(t) = C±B(t) and Y(t) = C 2 Z(t), where C 1 and C 2 are some arbitrary constants. We obtain 
another useful constraint on our functions from the (2, 2) component of the g equation of motion, 

-2 B ((ft + /3 2 X ) Z - B) + B 2 - B 2 (/3 0 + p x X - p) - (p 2 + p 3 X) Z 2 = ^ - 1. (44) 

It is easy to see that the r-dependent part of the left-hand side of Eq. (l44l) is s 2 (r). The right-hand side of Eq. (l44l) 
should also have the same dependence on r in order for the r dependence of both sides of the equation to cancel out. 
Therefore, in this case we have 


B' 2 s' 2 

^- 1 = ^- 1 = C 3 ^ 


(45) 


where C 3 is another arbitrary constant. From the (3,3) component of the equation of motion for we find 

/ / n 3 

s q — qs = 04 sq . 


(46) 


Equations (USD and (gH are the conditions which should be fulfilled by s and q for the consistency of the Einstein 
equations for both g and In summary, we find that the LTB- LTB combination is consistent only for particular 
subclasses of LTB metrics which satisfy all the conditions stated above. 


D. Bianchi I—FLRW combination 

We would also like to investigate the cases where one of the two metrics is homogeneous but anisotropic while 
the other one is both homogeneous and isotropic (i.e. has an FLRW form). Here we consider only Bianchi type 
I models, which are the simplest anisotropic models and capture most of the interesting anisotropic effects. The 
general properties of the cases where both metrics are anisotropic, simultaneously diagonal, and of the same Bianchi 
types within the Bianchi class A, which includes types I, II, VIo, VIIo, VIII, and IX, are discussed in Ref. [79}; we, 
therefore, do not consider those cases in this paper. As discussed in Sec. HI our main motivation for studying non- 
FLRW solutions in bigravity is the potential resolution of the problems with the standard scenario, in particular the 
instability issues, in this framework. However, there are also theoretical and observational arguments [851 ] in support 
of an anisotropic phase in the early Universe which approached isotropy at later times. It is, therefore, interesting 
also from this perspective to see whether such solutions are allowed in bigravity (see Ref. [79} for other motivations 
for studying anisotropies in bigravity, including a potentially interesting connection to dark matter). 

Bianchi type I, or simply Bianchi I, models are spatially homogenous and flat but the expansion rate is direction- 
dependent. In GR, these models have been studied for different sources with the equation of state p = cop. It has 
been shown in Ref. [13 that for cases with u> < 1 the anisotropic models evolve towards an FLRW universe, while for 
ui = 1 the process of isotropization does not take place. In the present study of Bianchi metrics we do not consider 
the question of isotropization and only investigate the consistency of such solutions in terms of the field equations 
and Bianchi constraint for both isotropic and anisotropic perfect-fluid sources. 

In this section we focus on the case where the physical metric is assumed to be of an anisotropic Bianchi I form 
while the reference metric is FLRW and flat; we call this case Bianchi I-FLRW. We show that for any choices of the 
matter source, isotropic or anisotropic, the Bianchi I -FLRW combination does not satisfy the conditions of the theory 


and is, therefore, not a consistent solution. 

The metrics for this particular case have the forms 

g^dx^dx" = — dt 2 + a 2 (t) dx 2 + a\ (t) dy 2 + a 2 (t) dz 2 , (47) 

f t _ lv dx tJ, dx 1 ' = —X 2 (t) dt 2 + b 2 (t) dx 2 , (48) 

where dx 2 = dx 2 + dy 2 + dz 2 , and x, y, and z are Cartesian coordinates. Here, ai, a 2 , and a 3 are the ^-metric 

scale factors along different directions, and b is the /-metric scale factor. The equation of motion for gives, 

independently of the matter source for the g metric, 

X 2 (Piaia 2 a 3 + p 2 (a\a 3 + aia 2 + a 2 a 3 ) b + /? 3 (ai + a 2 + a 3 ) b 2 + fab 3 ) - 3 bb 2 = 0, (49) 

V 2 (a 3 (a 2 (fa + faX) + b(fa+ faX)) + b (a 2 (fa + faX) +b(fa+ faX))) - X ( 2 bb + 6 2 ) + 2 bXb = 0, (50) 

V 2 (a 3 (ai (Pi + P 2 X) + b(p 2 + p 3 X)) + b (01 (P 2 + p 3 X) + b(p 3 + p 4 X))) - X (2 bb + b 2 ^j + 2 bXb = 0, (51) 

X 2 (a 2 (a 3 (Pi + P 2 X) + b (P 2 + p 3 X)) + b (ai (p 2 + P 3 X) + b (P 3 + p 4 X))) — X ^2 bb + b 2 ^j + 2 bXb = 0. (52) 
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From Eqs. and GD we see that ai = 02 , while Eqs. ED and El imply that 02 = 03 ; therefore ai = a 2 = a 3 . 
This means that the physical metric is forced to become isotropic. Therefore, using only the /-metric equation of 
motion and no other equations or conditions of the theory, we can conclude that one cannot have a Bianchi I physical 
metric while the reference metric maintains its FLRW form. Since / Mv is not sourced by matter, our conclusion is 
general and independent of what form the matter source takes. 


E. FLRW—Bianchi I combination 


1. Isotropic source 


Let us now assume that the reference metric takes the Bianchi-I anisotropic form while the physical metric is FLRW 
and flat; we call this case FLRW-Bianchi I. This means that each direction has a different scale factor in the / metric. 
The metrics, therefore, possess the forms 


g^dx^dx" = — dt 2 + a 2 ( t ) dx 2 , 

f^dx^dx” = —X 2 ( t) dt 2 + b\ ( t ) dx 2 + b 2 (f) dy 2 + b 2 ( t) dz 2 


(53) 

(54) 


where again dx 2 = dx 2 + dy 2 + dz 2 . Here, a is the scale factor for g^ v , and 61 , 62 and 63 are the /-metric scale factors 
along different directions. In addition, let us assume that the matter source, which couples to the physical metric g^, 
is a homogeneous and isotropic perfect fluid. The stress-energy-momentum tensor has, therefore, the form given in 
Eq. (1^1) . 

For the metric forms (l53l) and (l5Hl . the g- metric Einstein equation El reads 


3—T — P — Po + Pi 


bi 


- Pi 


b \ b 2 + 61/3 + &2^3 


P 3 


b\ 62 &3 


■P = , 


P= , 


■ Pi \X- 
Pi (x 


+ 2-+p = p 0 +pi(X + 


b 2 + 63 
a 

b\ + 63 

a 

bi + 6 2 
a 


■ Pi 

■ Pi 

■ Pi 


X (62 + 63) 62^3 


a 


a 2 

X (h + 63) 


hb 3 

a 


a 2 

X (hr + 62) 


6162 

a 


a 2 


+ P3 
+ p3 
+ P3 


Xbob, 


' 2 U 3 


a 

Xbib 


1«3 


a 

Xbib 


ivi 


cr 


(55) 

(56) 

(57) 

(58) 


It is straightforward to see from these equations that b±, 62 , and 63 are equal. This means that if we initially assume 
an anisotropic form for the reference metric while the physical metric and the matter source are both assumed to be 
isotropic, the structure of the Einstein equations automatically forces the reference metric to also be isotropic. 


2. Anisotropic source 


Let us now relax the isotropy condition on the matter source and study the FLRW- Bianchi I scenario when the 
physical metric is coupled to a source which is anisotropic and of a Bianchi I type. The stress-energy-momentum 
tensor in this case takes the form 


K 0 = -Pit), 

(59) 

£ 

II 

1 ? 

(60) 

T g1=P2(t), 

(61) 

T g3 =P3(t), 

(62) 


where the components of the fluid pressure, p±, P 2 , and P 2 , are allowed to be different along different directions and 
therefore create anisotropy in the fluid. 

It follows from the Bianchi constraint 0 that in this case the /-metric lapse X ( t ) is given by (here all the variables 
depend only on t ) 


1 pia 2 (61 + b 2 + b 3 )'+ p 2 a (bib 2 + frifr 3 + 62^3)'+ P3 (M2&3)' 
a 3 pia 2 + 2p 2 a (61 + b 2 + b 3 ) + p 3 (foib 2 + b 3 b 3 + b 2 b 3 ) 


( 63 ) 
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The Einstein equation for g^ gives 


O a _ a , a ft + ft + ft 

3—« — P — Po + Pi- 

a z a 




4 + 2- + Pi = A + A x + (A + AX) ?A- + (ft + AX )? 2 +1,3 


• 2 77 

^ + 2- + p 2 = ft + ftX + (ft + ftX) + (ft + ftX) 


• 2 7 T 

^ + 2- + p 3 = A) + ftX + (ft + ftX) + (ft + ftX) 


a 

b 3 + fti 
a 

62 + &i 
a 


(64) 

(65) 

( 66 ) 
(67) 


and the /-metric Einstein equation gives 

ft a 3 X 2 + fta 2 X 2 (61 + ft + ft) + ft«X 2 (61ft + ftft + ftft) + ftftftftX 2 — ftftft — 6162ft — ft ft ft = 0 , 

a 2 X 2 (Pi + ftX) + aX 2 (ft + ftX) (ft + ft) + 63 ^Xft + X 2 (ft + ftX) ft^ + ftXft — X ^ftft + ftft + ftft^ = 0 , 

a 2 X 2 (ft + ftX) + aX 2 (ft + ftX) (61 + ft) + 63 ^Xft + X 2 (ft + ftX) ft^j + ftXft — X ^ftft + 63 ft + ftft^ = 0, 

a“X 2 (ft + ftX) + aA’ 2 (ft + ftX) (ft + ft) + ft ^Xft + X 2 (ft + ftX) ft^ + ftXft — X ^ftft + ftft + ftft^ = 0. 

( 68 ) 


From Eqs. dMD -dSZD we see that because of the anisotropic matter source, ft, ft, and 63 satisfy different algebraic 
equations and can, therefore, be different from each other. Equations (IHSl) form a set of second-order differential 
equations with respect to ft, 62 , and 63 . Each equation has two independent solutions and, therefore, in spite of the 
fact that ft, 6 2 , and 63 satisfy the same differential equation, in general they can be different. We can, therefore, 
conclude that the combination FRLW-Bianchi I in the presence of an anisotropic matter source is in principle possible. 


IV. CONSISTENCY OF SOLUTIONS WITH ONLY ONE PERTURBED METRIC 

In the previous sections we explored various combinations of cosmologically interesting metric types for the two 
metrics of the theory of massive bigravity. Our investigation was, however, only at the level of the background 
dynamics of the Universe. In this section we study some other interesting cases where perturbations around the 
background metrics are considered. 

Perturbations are clearly crucial for cosmology and it is, therefore, interesting to ask whether there are any consis¬ 
tent, nonstandard ways to perturb the metrics. In all the previous works on the perturbative analysis of bigravity a 
standard recipe has been followed: the background metrics have been assumed to be FLRW and both the physical and 
reference metrics have then been perturbed around the FLRW solutions. As stated in Sec. HI our main motivation in 
this paper has been to try alternative solutions which might be free of various instabilities which appear in bigravity 
models at the level of perturbations. It would, therefore, be interesting if one could avoid the instabilities by finding 
alternative ways of perturbing the metrics. In addition, the usual perturbation equations are very complicated and it 
would be very helpful if one could find a way to simplify the equations. One of such approaches could be to work with 
solutions which leave the reference metric unperturbed while the physical metric is perturbed as usual. This scenario 
is physically justified because the reference metric does not couple to the matter components and cannot be measured 
directly from observations; the reference metric affects our observables only through its interactions with the physical 
metric. In what follows we analyze the consistency of this possibility using both scalar and tensor perturbations. 

Let us start with scalar perturbations up to linear order for both metrics. For simplicity we assume both metrics to 
be of the flat FLRW type. Including only scalar perturbations and using the notations of Ref. ( 68 j, the line elements 
for the perturbed physical and reference metrics g and have the forms 

ds 2 = a 2 (g) [— (1 + 24 < g ) eft 2 + 2diB g dx l dg + [(1 — 24> g ) + 2 didjE g \ eft*deft] , (69) 

ds 2 = b 2 (g) [— (1 + 24'/) X 2 cft 2 + 2diBfXdx l dg + [(1 — 24>/) ftj + 2ftftift] dx l dx^ 1 \ . (70) 

Here a and b are the scale factors corresponding to g and respectively, X is the lapse for and the 
perturbation quantities {4ft j, B g j, 4> g j, E g j} are allowed to depend on both conformal time g and space. Spatial 

indices are raised and lowered by the Kronecker delta. In this section a dot denotes a derivative with respect to the 

conformal time g. 
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In order to write down the expression for the perturbed matter stress-energy-momentum tensor, we assume a perfect 
fluid with an equation of state p = wp. In addition, we describe the matter perturbations with only one scalar field 
X- This procedure has been proposed in Ref. [87j] . With these assumptions and conventions we have 

6T 0 o = -(p + p)(3<f> g -E g>a -X,u), 

6T i 0 = -(p + p)x’ i , 

8T°i ={p + p)(B g ,i + x,i ), 



ST l j =w(p + p)(3$ g - E g p - x,uWj> 

(71) 

where we sum over indices l. Following Ref. [6Ej] and after some useful gauge fixing and transformations developed in 
Ref. [11|, we set $/ = x = 0 and arrive at the first-order perturbation equations 


2 n (3$g + fc 2 Eg) + a 2 ((1 + w)p{ 3$g + k 2 E g ) + rZ{ 3$ s + k 2 {E g - E f ))) 

+ 2 (fc 2 $g + - k 2 B g )) = 0, 

(72) 


2(X + 1 )$g + 2H(X + l)Tg - Zr{XB f - B g ) + (1 + w)p(l + X)B g = 0, 

2 (k 2 E g + 3$g) + 277(3Tg + 6$ 9 + 2 k 2 E g ) - 2 k 2 B g + 3 Za 2 rX(* f + tf ff ) 

(73) 


+ a 2 (-3(1 + w)p(2* g + iu(3$ g + k 2 E g )) + 2r(-3ZTg + Z{ 3$ s + k 2 (E g - E f )))) 

+ 2(9 U 2 - k 2 )^ g + 2fc 2 ($g - 2 HB g ) = 0, 

(74) 


E g — B g + 2 T~LE g + Za 2 r(E g — Ef ) — 4/g — 2 T~LB g + d^g = 0, 

(75) 


2 rU f k 2 E f -a 2 ZX 2 {k 2 E g - k 2 E f + 3$g) - 2rX'H f k 2 B f + 6 U 2 f rV f = 0, 

(76) 


2 U f T(X + 1 )T/ + Xa 2 Z(XBf - B g ) = 0, 

(77) 


rXE f -r{-2XUf + X)E f - X 2 (%r + rXVf + 2rH f B f + a 2 Z(E g - £») = 0, 

(78) 

where we have defined 



Z = Pi + 2 p 2 r + p 3 r 2 

(79) 


Z = Pi+ p 2 r (1 + X) + p 3 r 2 X , 

(80) 


- b oj h u - b 

r = -, U = ~ H f = 

a a b 

(81) 

Note that r in 
perturbed; i.e., 
find 

this section is not the radial coordinate. Let us now assume that only the physical metric g M! , is 
all the perturbative quantities for {T/,£?/, $/,£’/}, are vanishing. Looking at Eqs. ([75|)-([75|) we 


X 2 a 2 Z(k 2 E g + 3<Fg) = 0, 

(82) 


Xa 2 ZB g = 0, 

(83) 


X 2 a 2 ZE g = 0. 

(84) 


From Eqs. (I82p ~ (l84p it is clear that the g-metric scalar perturbations, B gi *I> 5 and E g , should vanish (which in turn 
implies, using Eq. m, that d'g should also be vanishing) unless the quantities Z and Z are both vanishing. In order 
to prove that the latter cannot be the case let us now assume Z = 0. In this case we have 

Pi + for (1 + X) + for 2 X = 0. (85) 

On the other hand, we know from the Bianchi constraint Q> that at the background level X satisfies the condition 

a b 


X = 

a b 

If we now insert the value of X from Eq. (1551) into Eq. (1551) we obtain 

. fo 


Pi ad + P 2 (ab + ba) + p 3 bb = 0 


fo 


2 (a 2 )' + p2 ( ab)'+ — (b 2 )' = 0. 


( 86 ) 


(87) 
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Finally, by integrating Eq. (1871) over conformal time we get 

foa 2 + 2foab +fob 2 = C => fo + 2 for + for 2 = (88) 

cr 

where C is an arbitrary constant. We, however, need to assume Z = 0 in this case, which then implies (7 = 0. The 
condition (1881) with (7 = 0 requires r to be a constant and given by a particular combination of the parameters /?i, /?2 
and /? 3 . From the background equations of bimetric gravity [68 ] one realizes immediately that this condition implies 
that the Universe is in a de Sitter state (or, trivially, that a = b = 0). 

Let us now study the opposite scenario, i.e. where only the reference metric is perturbed and the physical metric 
remains unperturbed. In this case the g-metric perturbations {'kg, B g} <f> 3 , E g } are vanishing. Now Eqs. (172l) - (175l) 
imply 


Zra 2 k 2 Ef = 0, 

(89) 

ZrXB f = 0, 

(90) 

3 Za 2 rX^! f - 2 Za 2 rk 2 E f = 0 

(91) 

Za 2 rEf = 0. 

(92) 


From Eqs. m-m we can again conclude that the /-metric scalar perturbations should also vanish in this case 
(note that >I >f is already vanishing due to our gauge choice), unless both Z and Z are vanishing. The latter case again 
corresponds to a de Sitter universe where r = const. We can, therefore, conclude, based on both cases studied here, 
that for any cosmological configurations different from a pure de Sitter universe, the possibility that only one of the 
two metrics is perturbed is excluded at the scalar level. 

A similar analysis can be done using tensor perturbations. The line elements for the physical and reference metrics 
in this case take the forms 


d 4 

ds) 


= a 2 ( 77 ) [-dr? + (5ij + ft?.) dx l dx j ] , 
= b' 2 ( 77 ) — X 2 dr] 2 + (i5 ij + h{^j dx l dx 9 


(93) 

(94) 


where ft?/ are the tensor perturbation quantities and are in general functions of space and conformal time 77. These 
quantities satisfy the relations 

h\ = 0, ft-* = 0, (95) 

for both ft 9 and ft/; here we again use the Kronecker delta to raise or lower spatial indices. 

In a perfect-fluid model there are no tensor perturbation modes in the stress-energy-momentum tensor, and as a 
result, the right-hand side of the perturbed Einstein equation for g^ vanishes. Because of the conditions (1^5|) we have 
only two degrees of freedom for each ft ;j - which are fully decoupled and satisfy the following first-order perturbation 
equations in Fourier space: 

ft's + 2 Uhs + k 2 h 9 + a 2 rZ (ft 9 - ft/) = 0, (96) 

ft> + ^2 H f - h f + k 2 X 2 h f + ^ Z (ft/ - ft 9 ) = 0. (97) 


As in the scalar case discussed above, we now assume that only the g metric is perturbed while the / metric is 
unperturbed (the argument is identical if we assume instead that only is perturbed); this means that ft/ = 0 . 
From Eq. (1!?T1) it is straightforward to see that the g-metric tensor perturbation ft 9 should also vanish, unless the 
quantity Z is vanishing. This corresponds to the condition (1881) , but contrary to the scalar case, here C is an arbitrary 
constant and is not required to vanish in order for the g and / tensor perturbations to decouple. One obvious solution 
is, however, again the case where (7 = 0, i.e. a purely de Sitter universe which does not correspond to the real 
universe. The only nontrivial class of solutions for which the tensor modes decouple while the scalar modes are still 
coupled is the case where (7 is nonvanishing. Whether or not these solutions are cosmologically interesting remains 
to be investigated; we leave this for future work@ 


2 The case for vector perturbations is not of interest to us here, but looking at the corresponding perturbative equations presented in 
Ref. fs7il it seems to us that the case where vector perturbations are nonvanishing only for the physical metric while the reference metric 
remains unperturbed is possible. 
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Metric combinations 

(physical metric) 

fp V (reference metric) 

7 n 9 

fMU 

Possibility 

Reference 

FLRW (k) 

FLRW (k) 

PF 

/ 

Standard 

FLRW (k g ) 

FLRW (k f ) 

PF 

X 

Present work and Ref. [55] 

FLRW 

Lemaitre 

PF 

X 

Present work 

FLRW 

Lemaitre 

Inhom. 

/ 

Present work and Ref. [26] 

Lemaitre 

FLRW 

Any 

X 

Present work 

FLRW 

LTB 

Any 

X 

Present work 

LTB 

FLRW 

Any 

X 

Present work 

LTB 

LTB 

PF 

/I 

Present work 

BianchiI 

FLRW 

Any 

X 

Present work 

FLRW 

BianchiI 

PF 

X 

Present work 

FLRW 

Bianchi I 

Aniso. 

/ 

Present work 

Bianchi Class A 

Bianchi Class A 

PF 

/ 

Refs. [30, 79] 

Perturbed FLRW 

Perturbed FLRW 

Perturbed PF 

/ 

Standard 

Perturbed FLRW (scalars) 

Unperturbed FLRW (scalars) 

Perturbed PF 

X 1 

Present work 

Unperturbed FLRW (scalars) 

Perturbed FLRW (scalars) 

Perturbed PF 

x b 

Present work 

Perturbed FLRW (tensors) 

Unperturbed FLRW (tensors) 

Perturbed PF 

XU 

Present work 

Unperturbed FLRW (tensors) 

Perturbed FLRW (tensors) 

Perturbed PF 

x c 

Present work 


a There are conditions which must be satisfied; see the text for details. 
b Unless the metrics are de Sitter; see the text for details. 

c Unless the metrics are de Sitter or satisfy a specific condition; see the text for details. 

Table I. List of all metric combinations studied in the present work, as well as a few other interesting combinations studied 
in the literature. / and x denote consistent and inconsistent cases, respectively. Here, “PF,” “Inhom.” and “Aniso.” stand, 
respectively, for isotropic and homogeneous (perfect fluid), isotropic but inhomogeneous, and homogeneous but anisotropic 
matter sources. “Any” stands for a matter source of any type, k is the spatial curvature of an FLRW metric; in cases with 
explicit indices g and / for k the two metrics are assumed to have different spatial curvatures. 


V. CONCLUSIONS 

One of the problems of contemporary cosmology is that there are very few, if any, well-studied, and well-motivated 
alternatives to the standard model that are viable and distinguishable from ACDM. Massive bigravity models 
particular model discussed in [H}) could be one of such rare cases if the problem of perturbation instability 
could be overcome. However, most studies of bigravity models confined themselves to homogeneous and isotropic 
metrics with no or identical spatial curvatures. Before ruling out bigravity cosmology one should, therefore, see if the 
problems could be solved or alleviated when a different geometry is chosen. This paper was devoted to addressing 
this question by scanning several relatively simple possible metric combinations in search for consistent cases. Future 
work will be necessary to actually solve the equations for some of such cases and identify cosmologic.ally viable ones. 

We, however, found that there are only a few metric combinations that survive our analysis. In most cases, two 
different metrics are impossible unless a suitable modification is made to the matter source beyond the standard 
perfect-fluid assumption. Our results are summarized in Table HI In particular we find that the only alternative 
combination with standard matter source is an LTB-LTB model (in addition to the Bianchi-Bianchi models previously 
studied in Refs. naizi), subject to some constraints. Combinations like FLRW-Lemaitre and FLRW-Bianchi I all 
require either inhomogeneous or anisotropic sources. We also investigated the question of whether having linear 
perturbations in just one metric is theoretically consistent. We found, not unexpectedly, that there are no consistent 
cases except for purely de Sitter backgrounds. 

Whether any of the surviving combinations give rise to viable cosmological models at both the background and 
perturbative levels is, however, entirely to be seen. 
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Appendix A: Impossibility of mapping Lemaitre to FLRW in bigravity 


In this appendix we prove that it is impossible to map the case where the physical metric v is FLRW and the 
reference metric has a lapse and scale factor depending on both r and t to a combination where both metrics are 
of an FRLW form. 

We start with the metrics in terms of the coordinates r and t, 


where 


g^dx^dx" = —dt 2 + a 2 (t) dx g 2 , 

(Al) 

f^dx^dx" = —X 2 (■ t , r) dt 2 + b 2 {t , r) dxf 2 , 

(A2) 

dx g 2 = —!——- + r 2 dtt 2 , 

1 — k g r z 

(A3) 

= — 2 + r 2 dVt 2 , 

1 — kfr z 

(A4) 


and dtt 2 = d, 6 2 + sin 2 (9) dcj) 2 . We now want to see if it is possible under any conditions to rewrite the metrics in an 
FLRW form when we transform the coordinates r and t to some new coordinates f and t. We assume the new and 
old coordinates to be related as 


t = T(t,f), (A5) 

r = R(t,f) . (A6) 

Under these transformations the g-metric scale factor a and the /-metric scale factor and lapse, b and X 1 should 
transform as 


a (t) —> a ( t , f) , (A7) 

X (t, r) X (f, r ) , (A8) 

b (t, r) -a 6 (t, f) . (A9) 

In order to know how the line elements (IA3I) and (lAdl) look under the coordinate transformations, we should first see 
how dt 2 and dr 2 transform. We have 


dt 2 = T 2 dP + 2 T'T dtdf + T , 2 df 2 , (A10) 

dr 2 = R 2 dP + 2 R'Rdtdr + R , 2 df 2 , (All) 


where an overdot denotes a derivative with respect to t and a prime denotes a derivative with respect to r. Using all 
these relations, the transformed line elements for the metrics read 


d~si = - T 2 - 


ds 2 f = - X 2 T 2 - 


1 - k g R 2 ' 
b 2 


R 2 ) dt 2 + 2 


1 -k f R 2 


R 2 I dt 2 + 2 


1 - k R R 2 
b 2 


R'R - T'T dtdf + 


1 - k g R 2 


R ' 2 - R ' 2 df 2 + a 2 R 2 dfl 2 , (A12) 


l-k f R 2 


R'R - X 2 T'T dtdf + 


b 2 


1 - k f R 2 


R - X T J df 2 + b z R 2 dn 2 . 

(A13) 


Now in order to have both metrics in an FLRW form, we need to set the following constraints on our transformed 
metric components. From Eq. (1A12I) we obtain 


r2 -i 

(A14) 

1 -k g R 2R,R - T ' T = 0 ' 

(A15) 

a 2 r ,2 rji /2 

1 - k g R 2 l _ k g f 2 ’ 

(A16) 

a 2 R 2 = B 2 (t)f 2 , 

(A17) 
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and from Eq. (IA13I) we find 


X 2 T 2 - 
b 2 

1 - k f R 2 


- _____ Tt 2 

1 -k f R 2 

R'R - X 2 T'T 


C 2 (t) , 

0 , 


n/2 _ y2rp/2 __ ^ ft) 

1-fc/i? 2 1-fc/f 2 ’ 

6 2 i? 2 = D 2 (t)r 2 , 

where A, B , C, and ZJ are arbitrary functions of only t. Now from Eqs. (IA17I) and (IA21I) . we get 

d 2 _ B 2 (t) 

b 2 “ D 2 ft) ’ 


(A18) 

(A19) 

(A20) 

(A21) 

(A22) 


and from Eqs. (1A15I) and (1A19I) . we get 


b 2 _ a 2 

(1 - k f R 2 )X 2 ~~ 1 - k g R 2 ' 

Now using the condition (1A23I) and Eqs. (IA14I) and (IA18I) . we find 


A 2 (i) 


c 2 (t) 

X 2 


(A23) 


(A24) 


This tells us that X must be a function only of t. Combining this with Eqs. (1A22I) and (I A23I) . we immediately see 
that R must also be a function only of f, and as a result R' = 0. Taking this into account, Eq. (1A19I) implies that 


T'T = 0. (A25) 

Let us first discuss the option where T = 0. In this case Eq. (IA14I) tells us that a must be a function only of t. Given 
that a and R are both functions only of t , Eq. (1A17I) immediately gives us a contradiction since the left-hand side is 
a function only of t, while the right-hand side is a function of both t and r. For the second option where we assume 
T' = 0, Eqs. (IA16I) and (IA20I) imply that D(t) = B(t) = 0, which of course means that we cannot have FR.LW metrics 
under any coordinate transformations. 
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